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Abstract 

Background: Reflexology is an alternative medical practice that produces beneficial effects by applying pressure to 
specific reflex areas. Our previous study suggested that reflexological stimulation induced cortical activation in 
somatosensory cortex corresponding to the stimulated reflex area; however, we could not rule out the possibility of 
a placebo effect resulting from instructions given during the experimental task. We used functional magnetic 
resonance imaging (fMRI) to investigate how reflexological stimulation of the reflex area is processed in the primary 
somatosensory cortex when correct and pseudo-information about the reflex area is provided. Furthermore, the 
laterality of activation to the reflexological stimulation was investigated. 

Methods: Thirty-two healthy Japanese volunteers participated. The experiment followed a double-blind design. Half 
of the subjects received correct information, that the base of the second toe was the eye reflex area, and pseudo- 
information, that the base of the third toe was the shoulder reflex area. The other half of the subjects received the 
opposite information. fMRI time series data were acquired during reflexological stimulation to both feet. The 
experimenter stimulated each reflex area in accordance with an auditory cue. The fMRI data were analyzed using a 
conventional two-stage approach. The hemodynamic responses produced by the stimulation of each reflex area 
were assessed using a general linear model on an intra-subject basis, and a two-way repeated-measures analysis of 
variance was performed on an intersubject basis to determine the effect of reflex area laterality and information 
accuracy. 

Results: Our results indicated that stimulation of the eye reflex area in either foot induced activity in the left middle 
postcentral gyrus, the area to which tactile sensation to the face projects, as well as in the postcentral gyrus 
contralateral foot representation area. This activity was not affected by pseudo information. The results also indicate 
that the relationship between the reflex area and the projection to the primary somatosensory cortex has a lateral 
pattern that differs from that of the actual somatotopical representation of the body. 
(Continued on next page) 
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Conclusion: These findings suggest that a robust relationship exists between neural processing of somatosensory 
percepts for reflexological stimulation and the tactile sensation of a specific reflex area. 

Keywords: Reflexology, Neuroimaging, Functional magnetic resonance imaging, Primary somatosensory cortex, 
Placebo effect 



Background 

Reflexology is an alternative medical practice that pro- 
duces beneficial effects on the human body by applying 
pressure to specific points or areas on the feet, hands, 
and ears called "reflex areas." It is a kind of massage 
therapy that is performed without using any special in- 
struments (i.e., uses hands only). The modern concept of 
reflexology was introduced in both the U.S. and Europe 
in the early 20th century [1], and it is believed to help 
eliminate stress, improve blood circulation, and restore 
the psychological balance of the body [2]. Although it is 
practiced throughout the world and several clinical stud- 
ies have shown positive effects of this practice, the pre- 
cise mechanisms underpinning its clinical effectiveness 
have not been elucidated [3-5]. Furthermore, the physio- 
logical basis of reflexology is not fully understood. Each 
reflex area is believed to correspond to an organ or part 
of the human body; however, this correspondence had 
not been sufficiently investigated. Recently, neuroscien- 
tific research has examined the scientific underpinnings 
of complementary and alternative medicine such as acu- 
puncture. Goldmann et al. [6] found that adenosine me- 
tabolism may influence the local antinociceptive effect of 
acupuncture, and Imai et al. [7] reported with acupunc- 
ture an improvement in autonomic function imbalance 
in rats under restraint stress. Previous neuroimaging 
studies using acupuncture have reported that acupoints 
related to different functions activate different areas of 
the primary somatosensory cortex (SI) [8] or different 
patterns of neural activity in the frontal-limbic- striatal 
regions [9]. Yoo et al. reported a difference in neural ac- 
tivity between real and sham acupunctural stimulation 
in the primary and secondary somatosensory cortex [10]. 
As for touch massage, including reflexology, a relaxing 
effect of reflexology was observed using electroencephal- 
ography (EEG) [11], pleasant touch massage induced ac- 
tivation of the pregenual anterior cingulate cortex [12], 
and reflexological treatments affected resting-state 
neural activity in the retrosplenial/posterior cingulated 
cortex [13]. 

In our previous study, we sought to provide experi- 
mental evidence of a somatotopical relationship between 
cortical activity in SI and sensory stimulation of reflex 
areas on the foot using functional magnetic resonance 
imaging (fMRI) [14]. We found that the area activated 
by stimulation of the reflex area was consistent with the 



somatotopic representation of the corresponding or a 
neighboring body part in SI. We interpreted these re- 
sults as evidence of a physiological relationship between 
the reflex area and a body part or organ. However, an al- 
ternative interpretation is that the findings were the re- 
sult of the placebo effect; it is possible that a sensory 
expectation resulting from a priori knowledge of the 
reflex area influenced the cognitive processing of the 
perceived reflexological stimulation. Blankenburg [15] 
reported that a somatosensory illusion called the "cuta- 
neous rabbit" [16] evoked cortical activity in SI that 
somatotopically corresponded to the illusory percepts. 
Furthermore, this illusion has been shown to occur out- 
side the body using a stick held across the left and right 
index fingers [17], indicating top-down cognitive pro- 
cessing such as might be observed in a body-object 
interaction when processing a tactile sensation. Thus, it 
is necessary to determine whether top-down modulation 
of somatosensory percepts resulting from a priori know- 
ledge of the reflex area influences the relationship be- 
tween the tactile sensation produced by stimulation of 
the reflex area and cortical activation in SI. 

In the present study, we used fMRI to investigate 
whether SI activation evoked by reflexological stimula- 
tion of a reflex area is affected by a priori information 
about the reflex area. We hypothesized that if stimula- 
tion of the reflex area activated SI representation of the 
body part corresponding to the reflex area regardless of 
the accuracy of the a priori information, then the 
somatotopical relationship between stimulation of the 
reflex area and neural activity in SI was robust. The eye 
reflex area was chosen as a stimulation point because it 
induced significant activation of the face representation 
area in the SI in our previous study [14], The eye reflex 
area is located in two places: at the base of the second 
and third toes. Thus, we were able to test the hypothesis 
in a reflex area corresponding to the same body part by 
assigning correct (eye reflex area) or pseudo- (shoulder 
reflex area) information to separate toes (Figure 1). 
Either the correct or the pseudo-information was given 
to participants just before the fMRI was performed. 

Additionally, we examined whether laterality of the 
foot reflex area was consist with the expected projection 
to SI. Our previous study [14] showed significant activa- 
tion of the face representation area in the left SI during 
stimulation of the eye reflex area on the right foot. 
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# All instructed areas are eye reflex area. 

# Half of the subjects were instructed 
the reverse assignment. 



Figure 1 Reflex area stimulated in the fMRI experiment. 

Detailed legend: Location of the eye reflex area on the base of the 
second and third toes and assignment of the correct- and pseudo- 
reflex-area information for each reflex area. 



However, the lateral relationship between stimulation of 
the reflex area on the contralateral foot and cortical acti- 
vation in SI has not been investigated. Thus, the present 
fMRI experiment had four experimental conditions: par- 
ticipants were stimulated at the two eye reflex areas on 
the right and left foot with each toe associated with ei- 
ther the correct- or pseudo-information condition. A 
two-factor factorial analysis of variance (ANOVA) was 
performed to clarify the effect of the accuracy of the 
assigned reflex- area information and the laterality of the 
cortical activation of the stimulated reflex area in SI. 

Methods 

Participants 

Thirty-two healthy Japanese volunteers participated in 
the present study (21 males and 11 females; mean age, 
22 ± 2 years, range, 20-31 years). All participants were 
judged to be right-handed according to the Edinburgh 
handedness inventory [18]. They had no history of 
neurological or psychiatric disease. It was confirmed that 
no participant had detailed knowledge of reflexology. All 
participants provided written informed consent, and the 
experimental protocol was approved by the Ethical 
Committee of the Tohoku University School of Medi- 
cine. The experiments were undertaken in compliance 



with national legislation and the Code of Ethical Principles 
for Medical Research Involving Human Subjects of the 
World Medical Association (Declaration of Helsinki). 

Experimental procedure 

The present experiment was a double-blind study. Be- 
fore the fMRI scan, each subject was acquainted with 
the reflex areas on the base of the second and third toe 
of their feet by an instructor. Half of the subjects were 
told that the base of their second toe was the eye reflex 
area and that the base of their third toe was the shoulder 
reflex area. Thus, those subjects received correct infor- 
mation about the reflex area for the base of the second 
toe and pseudo-information about the reflex area for the 
base of the third toe on both feet. The other half of the 
subjects received the opposite instructions, i.e., they re- 
ceived pseudo-information for the second toe and cor- 
rect information for the third toe on both feet. These 
instructions are summarized in Figure 1. After providing 
this information, the instructor marked each reflex area 
on the subjects feet. An experimenter and an assistant 
(who were not the instructor) entered the MRI scanner 
room together with the subject. The subject was placed 
in a supine position in the MRI scanner. To prevent 
head movements, the subjects head was fastened by a 
band, and the space between the head coil and the sub- 
ject s head was stuffed with sponge pads. A semi-lucent 
screen was put in front of the subjects face, and the sub- 
ject was instructed to gaze at the fixation cross on the 
screen throughout the fMRI measurement. The experi- 
menter wore a pair of headphones to receive an auditory 
cue, which was presented using Presentation software 
(Neurobehavioral Systems, Albany, CA, USA). The ex- 
perimenter stimulated each reflex area using a wooden 
stick with the right hand in accordance with the auditory 
cue. Because this was a double-blind experiment, the ex- 
perimenter and assistant did not know whether the cor- 
rect or the pseudo-information was associated with each 
reflex area. The strength of pressure stimulation was ad- 
justed according to each subjects self-assessment immedi- 
ately prior to the fMRI scan. The experimenter stimulated 
each reflex area on the subject in the MRI, and the stimu- 
lation strength was adjusted by ascertaining a comfortable 
level for each subject. The assistant held the subjects feet 
throughout the scanning procedure. 

The fMRI experiment was conducted using a block de- 
sign consisting of four sensory stimulation tasks that in- 
volved stimulating the base of the second or third toes 
of subjects' feet. The duration of each block was 5 s, 
with a 10-s interval between blocks. The experimenter 
continuously stimulated one reflex area in each task 
block. Each task was repeated 10 times; a 26-s rest 
period was imposed before the first block, and a 24-s 
rest period was imposed after the last block. Thus, the 
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total time of fMRI scanning was 10 min 50 s. The order 
of tasks was counterbalanced among the subjects. 

fMRI data acquisition 

All images were acquired using a Philips Achieva 3T 
MRI scanner (Philips, Amsterdam, The Netherlands). 
The fMRI time series data covering the entire brain was 
acquired using a T2-weighted gradient echo echo-planar 
imaging (GE-EPI). The parameters of the experiment 
were as follows: repetition time (TR), 2,500 ms; echo 
time (TE), 30 ms; flip angle (FA), 80° 42 slices; field of 
view (FoV), 192 x 192 mm; 64 x 64 matrix; slice thick- 
ness, 3 mm; slice gap, 0 mm. We acquired 260 scans per 
subject. To acquire a fine structural whole-head image, 
magnetization-prepared rapid-acquisition gradient-echo 
(MP-RAGE) images were obtained (TR, 6.5 ms; TE, 
3 ms; FA, 8° FoV, 240 x 240 mm; 240 x 240 matrix; 162 
slices; slice thickness, 1.0 mm). 

Data analysis 

Data preprocessing and statistical analyses of the fMRI 
data were carried out using statistical parametric map- 
ping software (SPM8; Wellcome Trust Center for 
Neuroimaging, London, UK). The effect of head motion 
across the scans was corrected by realigning all images 
to the first one. No data were excluded from subsequent 
analyses based on our head-motion exclusion criteria 
(i.e., head motion in the fMRI scanner greater than 3 
mm). The structural image volume was then co- 
registered with the mean image of realigned echo-planar 
image (EPI). All EPI were spatially normalized to the 
Montreal Neurological Institute (MNI)-Tl template 
using the parameter to co-register and to normalize the 
structural image for the MNI-T1 template obtained by a 
segmentation process for each subject. Finally, each scan 
was smoothed with a Gaussian filter in a spatial domain 
(8-mm full-width at half maximum). 

The fMRI data were analyzed using a two-stage 
approach in SPM8. In the first-level analysis, the 
hemodynamic responses produced by different experi- 
mental conditions were assessed at each voxel using a 
general linear model on an intrasubject basis. The ana- 
lysis was based on the hypothesis that the hemodynamic 
responses were induced by the stimulation of the eye 
reflex area on the right foot with correct information 
[R(correct)] and pseudo [R(pseudo)] and the stimulation 
of the eye reflex area on the left foot with correct 
information [L(correct)] and with pseudo-information 
[L(pseudo)]. The hemodynamic response was assumed 
to be the canonical hemodynamic response function of 
5-s duration. Global intensity normalization was 
performed on the average within-brain fMRI signal in 
each image, and low-frequency confounding effects were 
removed using a high-pass filter with a 128-s cutoff. A 



multiple regression analysis was performed on each 
voxel to detect the regions where MR signal changes 
were correlated with the hypothesized model to obtain 
the partial regression coefficients of each voxel. 

The second-level analysis was performed using a two- 
way repeated-measures ANOVA on an intersubject 
basis. One factor was the laterality of the stimulated re- 
flex area (left or right; two conditions), and the other 
factor was the accuracy of the assigned reflex-area infor- 
mation (correct or pseudo; two conditions). Further- 
more, we performed a subtraction analysis of the main 
effect, which showed significant differences and both 
common and respective activation in the four experi- 
mental conditions. To depict an activation peak in the SI 
region, an anatomical mask image of the bilateral 
postcentral gyrus was used as an inclusive mask in the re- 
sults section of SPM together with the whole-brain ana- 
lysis. The mask image was generated using the SPM 
Anatomy toolbox [19]. The statistical threshold was set at 
p < 0.05 (corrected for family-wise error (FWE) by voxel 
level). Additionally, to assess the potential influence of sex 
differences, a three-way repeated-measures ANOVA with 
the factors stimulated reflex area laterality, accuracy of the 
assigned reflex-area information, and sex (male or female; 
two conditions) was performed. The main effect of sex 
and the interaction with other factors were tested with a 
statistical threshold of p < 0.05, corrected for FWE. 

To examine activation induced under each experimen- 
tal condition on an individual level, we created plots of 
the average activity levels within the regions of interest 
(ROIs) across subjects to evaluate local signal changes in 
the activated regions. This analysis was carried out using 
the MarsBaR toolbox [20]. ROIs were 4-mm-radius 
spheres centered within activation foci obtained by the 
second-level analysis. Because we focused on SI activity, 
an anatomical mask image of the bilateral postcentral 
gyrus was used to extract activation foci from the statis- 
tical analysis. The mask image was the same image used 
in the second-level analysis. Additionally, a portion of 
shoulder representation in SI [x = 34, y = -48, z = 70] 
that showed significant activation when the shoulder re- 
flex area was stimulated in our previous study [14] and 
its contralateral side [x = -34, y = -48, z = 70] were used 
as ROIs to determine whether the pseudo-information 
stimulated the shoulder reflex area. A time series of per- 
cent signal changes in each ROI per scan from the be- 
ginning of each experimental condition (5 s) to the next 
resting condition (10 s) was extracted from preprocessed 
fMRI data on each subject, and the average time series 
of percent signal change across subjects was calculated. 

Results 

The results of the second-level two-way repeated mea- 
sures ANOVA are summarized in Table 1. The main 
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Table 1 Cortical areas showing significant differences in the two-way repeated-measures ANOVA 

Area Brodomann's area MNI coordinate F-score 

x y z 

Main effect of truth of assigned reflex area information 



N.S. 

Main effect of laterality of reflex area 



L. middle cingulate gyrus 


BA6 


-8 


-4 


46 


61.69 


R. supplementary motor area 


BA6 


8 


6 


48 


60.58 


L. paracentral lobule 


BA4a 


-6 


-38 


66 


374.94 




BA3b 


-8 


-40 


74 


184.84 


R. paracentral lobule 


BA4a 


6 


-28 


72 


276.24 




BA3b 


10 


-44 


70 


173.63 


L. insula 




-34 


-24 


18 


42.42 


R. insula 




32 


-24 


18 


37.39 


L. Cerebellum (lobule V) 




-14 


-40 


-24 


125.91 


R. Cerebellum (lobule V) 




14 


-42 


-24 


156.89 



Interaction between two factors 
N.S. 

The activation peak in bold was observed within SI area defined by the Anatomy toolbox [20]. 



effect of reflex-area laterality showed a significant 
bilateral difference in a broad region that included the 
superior SI, where the somatotopic representation of the 
foot is generally found, and the adjacent regions of the 
precentral gyrus (BA 4a) and superior parietal lobule. 
Additionally, the bilateral medial prefrontal region, in- 
sula, and cerebellum also showed significant differences. 
In contrast, no significant difference was found for the 
main effect of correct/pseudo-information or the inter- 
action between information and laterality. Figure 2 
shows the cortical regions where a significant difference 
in the main effect of laterality was observed in superior 
SI and an average time series of the percent signal 
changes of subjects for each SI activation peak at the on- 
set of each experimental condition. The cortical areas 
significantly activated in each subtraction analysis are 
shown in Table 2. From the direct comparison of the 
left-foot versus the right-foot stimulation, significant 
activations were observed in the superior part of SI and 
in the neighboring motor or parietal region and insula 
in the stimulated foots contralateral hemisphere and 
in the ipsilateral cerebellar anterior lobe (Lobule V). 
Conversely, subtraction analysis between reflexological 
stimulation and true versus pseudo-information did not 
reveal any significant activation. The bilateral portions of 
the shoulder representation in SI were not significantly 
activated by stimulation with pseudo-information even 
though the statistical threshold was set at a more liberal 
threshold (p < 0.05, uncorrected). 

The cortical areas that showed significant common ac- 
tivation across the four conditions are summarized in 



Table 3, and the activation peaks within the SI region 
obtained by each of the four experimental conditions are 
summarized in Table 4. Reflexological stimulation of 
each foot reflex area activated the superior part of the 
contralateral SI hemisphere. Stimulation of both eye re- 
flex areas significantly activated the left middle part of 
the postcentral gyrus, where the somatotopic representa- 
tion of the face is generally found. However, the homolo- 
gous contralateral region of the middle part of the 
postcentral gyrus was not significantly activated by 
stimulation of the eye reflex areas even though the 
statistical threshold was set at more liberal threshold 
(p < 0.05, uncorrected). Moreover, we observed signifi- 
cant bilateral activation of the paracentral lobule, located 
slightly lateral to the area where the main effect of reflex 
area laterality was located, and observed that the activa- 
tion focus was generally the somatotopic representation 
of the lower limb. Moreover, significant activation in the 
bilateral motor region, secondary somatosensory area 
(SII) within the rolandic operculum, insula, and thal- 
amus was observed. Figure 3 shows the common activa- 
tion peak in the SI area observed under all conditions 
and an average time series of the percent signal changes 
in each activated area and additional ROIs: the contralat- 
eral homologous region of the left middle part of the 
postcentral gyrus and the bilateral shoulder representa- 
tion in SI determined by the location that showed 
increased activity in response to stimulation of the 
shoulder reflex area in our previous study [14]. A 
stimulation-related signal change was also observed in 
the averaged individual time series data in the bilateral 
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-L(True) 
— L(False) 
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— L(False) 
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— R(False) 
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[ 10, -40, 76] 

Figure 2 Cortical areas showing a significant main effect of reflex-area laterality in SI and its activation plot. Detailed legend: Results of 
second-level analysis showing a significant main effect of reflex area laterality [F-contrast of {L(correct) + L(pseudo)} versus {R(correct) + R(pseudo)}] 
superimposed on axial sections of the MNI single-subject template. The crosshairs show the activation focus peak position on the (A) left and (B) 
right paracentral lobules. Each plot shows an average time series of percent signal change across the subjects on the (A) left and (B) right 
paracentral lobules. The horizontal axis represents time elapsed from the onset of each condition. The error bars show the standard error of 
the mean. 



Table 2 Cortical areas showing significant activation in each subtraction analysis 



Area 



Brodomann's area 



MNI coordinate 

y 



T-score 



Stimulation of reflex area on left foot: {L(correct) + L(pseudo)} - {R(correct) + R(pseudo)} 



R. supplementary motor area 
R. paracentral lobule 

R. rolandic operculum 
R. insula 
R. thalamus 

L. cerebellum (lobule V) 



BA6 

BA4/6a 

BA3b/4a 



8 

44 
32 
16 
-14 



Stimulation of reflex area on right foot: {R(correct) + R(pseudo)} - {L(correct) + L(pseudo)} 
L. middle cingulate gyrus BA6 -8 

L. paracentral lobule BA4a/6 -6 

BA4 -6 
BA3b/4a -10 
L. insula -34 
R. cerebellum (lobule V) 14 
Stimulation of reflex area with true information: {L(correct) + R(correct)} - {L(pseudo) + R(pseudo)} 
N.S. 

Stimulation of reflex area with pseudo information: {L(pseudo) + R(pseudo)} - {L(correct) + R(correct)} 
N.S. 



-28 
-40 

-32 
-24 
-24 
-40 

-4 
-38 
-26 
-44 

-24 
-42 



48 
72 
74 

22 
18 
0 
-24 

46 
66 
70 
70 

18 
-24 



7.78 

16.62 

15.32 

5.10 
6.11 
5.08 
11.22 

7.85 
19.36 
15.37 
14.07 

6.51 
12.53 



The activation peak in bold was observed within SI area defined by the Anatomy toolbox [20]. 
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Table 3 Cortical areas showing common activation under four experimental conditions 



Area 


Brodomann's area 


X 


MNI coordinate 

y 


z 


T-score 


Common activation to four conditions: L(correct) + L(pseudo) + R(correct) + R(pseudo) 








L middle cingulate gurus 




-14 


-32 


40 


5.38 


L posterior cingulate gurus 




-8 


-44 


20 


5.98 


L. postcentral gyrus 


BA3b 


-14 


-44 


64 


7.60 




BA1 


-60 


-22 


44 


5.09 


R. postcentral gyrus 


BA3b 


14 


-44 


64 


6.76 


L. rolandic operculum 




-48 


-22 


16 


8.53 






-54 


-2 


8 


7.99 


R. rolandic operculum 




54 


0 


8 


7.19 






48 


-26 


20 


7.02 


L. insula 




-38 


0 


12 


9.39 






-34 


-20 


16 


9.57 


R. insula 




38 


8 


12 


9.59 






34 


-22 


14 


9.57 


L. thalamus 




-8 


-10 


6 


7.67 


R. thalamus 




18 


-12 


10 


8.39 


L. Amygdala 




-20 


-2 


-6 


5.33 


R. Amygdala 




20 


-4 


-6 


5.10 



The activation peak in bold was observed within SI area defined by the Anatomy toolbox [20]. 



Table 4 Comparison of SI activation peaks observed 
under four experimental conditions 



Area Condition 


X 


MNI coordinate 

y 


z 


T-score 


Left superior part of SI area 










L(Jrue) 


-18 


-42 


60 


5.21 


L(pseudo) 


-18 


-42 


60 


5.49 


R(True) 


-12 


-48 


71 


9.88 


R(Pseudo) 


-12 


-44 


66 


9.50 


Right superior part of SI area 










L(True) 


14 


-42 


70 


8.03 


L(pseudo) 


12 


-44 


68 


8.22 


Rdrue) 


16 


-44 


60 


4.98 


R(Pseudo) 


16 


-40 


65 


5.24 


Left middle part of SI area 










L(True) 


-60 


-22 


44 


4.27 


L(pseudo) 


-58 


-24 


48 


4.87 


Rdrue) 


-58 


-22 


44 


4.82 


R(Pseudo) 


-60 


-26 


49 


4.96 



The SI area mask image was created using the Anatomy toolbox [20]. The 
statistical threshold was set at p < 0.05, corrected for FWE. The activation peak 
in italics was observed under a liberal threshold (p < 0.001, uncorrected). 



paracentral lobule and left middle SI, whereas the right 
contralateral portion of the middle SI and expected SI 
shoulder representation did not show any significant sig- 
nal changes. 

The results of the second-level three-way repeated- 
measures ANOVA, which assessed the influence of sex 
differences, revealed no significant differences for the 
main effect of sex, in the two-way interaction between 
sex and correct/pseudo-information and between sex 
and laterality, or in the three-way interaction among sex, 
correct/pseudo-information, and laterality. 

Discussion 

The present study aimed to test whether SI neural activ- 
ity induced by tactile stimulation at reflexological stimu- 
lation sites was affected by a priori knowledge about 
reflexology. The results suggest that a robust relation- 
ship exists between tactile stimulation of the eye reflex 
area and neural activity representing the somatosensory 
percept of the reflexological stimulation of the corre- 
sponding body part. A second-level ANOVA revealed no 
significant effect of pseudo-information about the reflex 
area. The middle part of the postcentral gyrus, which is 
associated with face representation in SI [21-23], showed 
significant common activation by all of the eye reflex 
areas. Thus, neural activity corresponding to the percep- 
tion of reflexological stimulation was not affected by 
a priori knowledge about the reflex area. These results 
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Figure 3 Cortical areas showing common activation foci under the four conditions and activation plots. Detailed legend: Results of 
second-level analysis showing common activation foci under the four conditions [L(correct) + L(pseudo) + R(correct) + R(pseudo)] superimposed 
on axial sections of the MNI single-subject template. Crosshairs show the activation focus peak position on the (A) left and (B) right paracentral 
lobules and the (C) left middle part of the postcentral gyrus. Crosshairs at (D), (E), and (F) show the center point of each designated ROI. The 
statistical threshold was set at p < 0.05, corrected for family-wise error. Plots show an average time series of percent signal changes across 
subjects on the (A) left and (B) right paracentral lobules and the (C) left middle part of the postcentral gyrus. Plot (D) shows an average time 
series of percent signal change across subjects on the contralateral coordination (C). Plots (E) and (F) indicate each average time series of percent 
signal changes across subjects on the bilateral shoulder representation in SI determined by the location obtained from stimulation of the 
shoulder reflex area in our previous study [14] and its contralateral coordination. The horizontal axis represents time elapsed from the onset of 
each condition. The error bars show the standard error of the mean. 



support our previous finding that reflexological stimula- 
tion induces a somatosensory process corresponding to 
the stimulated reflex area [14]. Several neuroimaging 
studies of acupuncture have also suggested that stimula- 
tion of different acupoints on upper and lower limbs 
led to activation in somatotopically different locations 
[8,24]. However, that relationship basically showed a 
connection between the acupoint location on a body 
part and the SI somatotopic representation; however, an 
acupoint does not necessarily have a clear relationship 
with a specific body part. Laser stimulation of a specific 
acupoint empirically related to ophthalmic disorders re- 
portedly induces visual association cortex activation 
[25]. Therefore, these differential findings may imply 
that different neural mechanisms are involved in pro- 
cessing perceptions of reflexological and acupunctural 
stimulation. 



Furthermore, our results indicate that the reflexolo- 
gical stimulation of the eye reflex area induces 
left-lateralized activation in the middle part of the 
postcentral gyrus, as the contralateral side (right) was 
not significantly activated. Anatomically, the final projec- 
tion of face tactile sensation is to the middle part of the 
postcentral gyrus of the contralateral hemisphere via the 
trigeminal nerve. Neuroimaging studies have shown 
functional lateralization of face representation in SI. For 
example, Eickhoff et al. [21] reported that unilateral 
stimulation of the face induced bilateral activation of the 
middle part of SI. Another fMRI study demonstrated 
that the lateral relationship between tactile stimulation 
and SI activity differed depending on which part of the 
face is stimulated [22], and Moulton et al. [23] suggested 
that face representation in SI is arranged segmentally by 
subdivision of the trigeminal nerve. These discrepancies 
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may be explained by the present results, which indicate 
that the somatotopic relationship between the reflexolo- 
gical stimulation of the eye reflex area and the neural 
processing of tactile sensations emanating from the cor- 
responding body part is characterized by a lateral pattern 
that differs from that of the actual relationship and that 
the stimulation of the eye reflex area was perceived as a 
tactile sensation emanating from a more specific part of 
the face. Two possible interpretations of the left- 
lateralized SI activation can be offered. First, it may be 
that all of the eye reflex areas we stimulated correspond 
to a unilateral face representation and that the reflex 
area corresponding to the contralateral eye is expressed 
in another area. Alternatively, tactile sensations activated 
by reflexological stimulation of the eye reflex area may 
be processed primarily in the left hemisphere of SI. It is 
necessary to examine the relationship between SI activity 
and the stimulation of reflex areas corresponding to 
other body parts to determine whether either of these 
interpretations is valid. Balslev et al. [26] reported that 
cortical activation-related eye muscle proprioception was 
located at a slightly anterior portion of the middle 
sensorimotor cortex area compared with our results 
(x, y, z = -36, -16, 40 for right eye and x, y, z = -38, -12, 
42 for left eye, on left sensorimotor cortex) and that the 
activation peak specifically for eye movement was lo- 
cated in a more anterior location. The present study did 
not investigate the relationship between reflexological 
stimulation and induced muscle activation on the corre- 
sponding body part because we did not measure related 
muscle activation, such as eye movement, using electro- 
myography or eye tracking. However, specific cortical 
activation-related eye movement, such as in the frontal 
eye field, was not observed in the present results. In con- 
trast with the results reported by Balslev et al. [26], we 
observed that activation of the middle SI might be 
reflected in activity induced not by voluntary eye move- 
ment but by the perception of stimulation. 

Activation of another SI location that corresponded to 
the stimulation of each foot was observed. Stimulation 
of the eye reflex area on each foot induced significant 
activation in the contralateral hemisphere of the primary 
sensorimotor cortex, medial prefrontal region, SII, in- 
sula, and the ipsilateral hemisphere of the anterior part 
of the cerebellum (lobule V). This relationship between 
somatosensory stimulation and laterality was already 
suggested in a previous study of somatosensory process- 
ing [27]. Tactile sensations of the foot project to the 
superior part of SI [28], and our previous study demon- 
strated a consistent relationship between reflexological 
stimulation of the left foot and cortical activation in the 
right superior part of SI [14]. Furthermore, the activated 
region spread to neighboring areas, including the pri- 
mary motor area (area 4a/6) and superior parietal lobule. 



Yoo et al. [10] reported that real acupunctural stimula- 
tion induced broader activation in the sensorimotor 
cortex and posterior parietal region and area 3a was se- 
lectively activated by real acupunctural stimulation as 
compared with sham acupunctural stimulation. Al- 
though the present result reflected a similar trend in the 
activation pattern, only a few overlapping areas were ob- 
served on a probabilistic map of area 3a, as defined by 
the Anatomy toolbox. This difference might be due to a 
difference in the stimulation procedure used by reflex- 
ology and acupuncture. As area 3a is responsive to 
proprioceptive sensation and reportedly responds to 
kinesthetic and motor stimulation [29], foot kinesthetic 
responses induced by reflexological stimulation might 
differ from those induced by acupunctural stimulation. 
The activation foci located in a more lateral portion of 
the superior SI were observed during bilateral stimula- 
tion of both feet, although the activation plot showed 
that the activity on the contralateral side tended to be 
greater than that on the ipsilateral side. In terms of the 
somatotopic representation of the foot in SI, the activa- 
tion focus becomes more lateral when the more anterior 
portion of the leg is stimulated [30,31]. In our experi- 
ment, the assistant held the subjects' feet in the same 
position during fMRI measurement. When the reflex 
area was stimulated, the assistant should have kept the 
subjects' feet in the same position as in the resting state; 
indeed, the activity in this area reflected the tactile sen- 
sation that arose from foot-position fixation. 

Our results showed a robust somatotopic relationship 
between stimulation of the eye reflex area and activation 
of the middle part of SI, with no effect of a priori infor- 
mation regarding the reflex area. However, our study 
had several limitations. We were unable to specify the 
exact relationship between a tactile sensation of the re- 
flex area and a specific part of the face. We were also 
unable to elucidate the mechanism underlying the bene- 
ficial effect of reflexological treatment because our 
experimental setup was focused on the relationship be- 
tween SI cortical activation and erroneous reflexology 
information. Further, the present study did not clarify 
the relationship between reflexological stimulation and 
induced muscle activation on eye or shoulder movement 
because we did not measure muscle activation using 
electromyography or eye tracking. As for the integrated 
study of neuroscience and acupuncture, we clarified a 
metabolic contributor to the local antinociceptive effect 
of acupuncture [6]. Reflexological stimulation induces 
changes in resting-state neural activity [13], and it is 
possible that this kind of change in resting-state neural 
activity might be a potential confounder of our experi- 
mental results. Thus, a multifaceted, systematic investi- 
gation that includes neuroimaging approaches is needed 
to reveal the physiological basis of reflexology. 
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Conclusions 

The relationship between cortical activity and reflexological 
stimulation is robust, even when a priori pseudo- 
information about the reflex area is given to subjects. The 
results of the present study support our previous finding 
showing that reflexological stimulation of the eye reflex 
area induced significant activity in the face representation 
area in SI of the left hemisphere and in the foot representa- 
tion area in SI of the contralateral hemisphere. 
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